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Abstract 

We consider the role of supersymmetric flat directions in reheating the Universe after 
inflation. One or more flat directions can develop large vevs during inflation, which can po- 
tentially affect reheating by slowing down scattering processes among inflaton decay products 
or by coming to dominate the energy density of the Universe. Both effects occur only if flat 
directions are sufficiently long-lived. The computation of their perturbative decay rate, and a 
simple estimate of their nonperturbative decay have led to the conclusion that this is indeed 
the case. In contrast, we show that flat directions can decay quickly through nonperturbative 
channels in realistic models. The mass matrix for MSSM excitations around flat directions 
has nondiagonal entries, which vary with the phase of the (complex) flat directions. The 
quasi-periodic motion of the flat directions results in a strong parametric resonance, leading 
to the rapid depletion of the flat direction within its first few rotations. This may preclude 
any significant role for the flat directions in reheating the Universe after inflation in models 
in which the inflaton decays perturbatively. 
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1 Introduction 



One of the most attractive mechanisms for producing the baryon asymmetry of the Universe, 
is the out-of-equilibrium decay of coherent scalar field oscillations along nearly F- and D-flat 
directions of the scalar potential in supersymmetric theories [U [2] . The minimal supersym- 
metric standard model contains many of these flat directions [3] and one or more can be 
expected to be excited during inflation [U [5] . There are many consequences of the evolution 
of large vevs along flat directions [6]. If, in the context of a grand unified theory (GUT), there 
exists a scalar operator O, which violates baryon number and for which (O) ^ 0, the CP 
violating decay of the flat direction will produce a large baryon asymmetry. Alternatively, 
any operator violating lepton number could be used to generate a net lepton asymmetry [7] 
which is subsequently converted to a baryon asymmetry by sphaleron processes. 

Neglecting nonperturbative effects, one typically finds that the scalar field oscillations 
of a flat direction persist after the decay of the inflaton [5j [8j [9]. Indeed, quite generally, 
the Universe becomes dominated by flat direction oscillations until these decay. It has also 
been recently argued [9] that flat directions can delay the thermalization of the inflaton 
decay products (resulting in a low reheat temperature) by providing a large mass to gauge 
bosons through their vev, and thus suppressing the rates of the scattering processes needed 
for thermalization. 

Here, we examine the effects of preheating on the fate of flat directions. We find that 
nonperturbative effects can be important and lead to the rapid decay of the flat directions, 
long before they can dominate the expansion rate of the Universe. Similarly, the rapid decay 
may preclude any delay in the thermalization process after the decay of the inflaton. These 
effects are duetto mixing between different MSSM excitations around many different sets of 
flat direction 



a 



and have been overlooked in the literature so far. 
We will work in the context where the inflaton, x/j, couples to standard model fields only 
through gravity 3 We will also assume for simplicity that there is only a single mass scale 
associated with the inflaton potential set by the COBE normalization [11] , with ~ Hi ~ 
10~ 7 Mp, where Hi is the Hubble parameter during inflation and Mp is the Planck scale. 
We also assume that after inflation, oscillates with a Planck scale amplitude. The energy 
density in oscillations can be written as 

Pi> = mlij 2 = mlM 2 P {R^/Rf (1) 



^wo or more mutually non-exclusive flat directions are required. As we show below, only in exceptional 
cases does a single flat direction exclude all others from being excited. 

2 Reheating in the context of nonperturbative inflaton decay has been discussed in |10j . 
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where is the Robertson- Walker scale factor when oscillations begin. For R > R^, the Uni- 
verse expands as a matter dominated Universe with Hubble parameter, H ~ m^(R^/R) 3 ^ 2 . 
Inflatons decay at Rd^p — {Mp/m^Y'^R^ when their decay rate 1^ = rn^/Mp ~ H. After 
inflaton decay, the Universe is radiation dominated and p r ^ ~ m 2 / 3 Mp°^ 3 (R^/R)' 1 , where 
p r ^ is the energy density in the radiation produced by inflaton decay 

Thermalization occurs when the relevant scattering rates become comparable to the Hub- 
ble parameter. Naive estimates indicate that this may be delayed due to suppressed cross 
sections of the form a ~ (a 2 /m 2 )(i?/i? rf ^) 2 where a is a characteristic coupling strength 
(presumably of order the gauge fine structure constant) and the suppression is due to the 
initially large particle energies produced in the decay (E ~ m^) which redshift to lower 
energies as the Universe expands. The resulting reheat temperature in this case would be 
expected to be low, Tp ~ 10 4 GeV [5]. However it was shown [12] that forward scattering in 2 
— > 3 interactions, though of higher order in a, lead to rapid thermalization. It was estimated 
that the inelastic rates were roughly Ti ne i ~ a(Mpjm^)n^a in which case thermalization 
occurs almost immediately after decay and restores the simple estimate Tp ~ (T^Mp) ~ 10 8 
GeV. This will be true so long as a 3 > m^/Mp. 

In supersymmetric theories, quantum fluctuations can drive several scalar fields to large 
expectation values along flat directions. In fact, in supergravity models with a Heisenberg 
symmetry [13], such as no-scale supergravity [H], we expect to generate vev's of order 
0o ~ Mp for any flat direction not involving stops [T5] . This vev may be smaller if the flat 
direction is regulated by non-renormalizable operators [8]. 

When the Hubble parameter after inflation is of order the supersymmetry breaking scale 
m^,, sfermion oscillations along the flat direction will begin. This occurs at R — R<f> ~ 
{m^lm^) 2 l 3 R^. The energy density in the oscillations of the flat direction field characterized 
by cf) is given by 

/V = m 2 ^ 2 = ml^KR^/R) 3 ~ m\ <f>l{R^/Rf (2) 

Typically, oscillations of <p will begin before inflaton decay, since (Rd^/R^) ~ Mp 3 m 2 J 3 /m 2 ~ 
10 3 . 

We note that although most standard model (SM) fields are massive when the flat di- 
rection is excited, inflaton decay is unimpeded if the inflaton can decay directly into the 
light fields making up the flat direction. Nevertheless, the scattering rates will be suppressed 
due to the large masses of exchanged particles [9]. To compute the thermalization rate of 
the decay products, we once again consider the inelastic 2^3 scatterings by replacing a 
with (a 2 /0 2 ) and 2 = (0gm 2 /m 2 i )(i?^/i?) 3 . In this case, thermalization is not suppressed 
unless a < ml 3 jm! 3 M P . For a 2 ~ 1CT 3 , one finds again a reheat temperature of order 
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Tr ~ a^^m^Mp/m^ ~ 10 8 GeV. Note that we have assumed 0o — Mp. For lower 0o, 
reheating at this stage will always be given by the instantaneous rate yielding Tp ~ 10 s 
GeV. However as we will see, the energy density of the Universe will become dominated by 
the oscillations alongthe flat direction and final thermalization will depend on the decay of 
the flat direction veMJ. 

If inflaton decay occurs before the decay of flat direction, the radiation produced by in- 
flaton decay will redshift faster than the flat direction oscillations. It is of interest, therefore, 
to compare the energy density in inflaton decay products with that in the flat direction 
oscillations at the time of decay (at = 0^). The ratio of energy density in radiation to 
oscillations is 



The oscillations end when the decay rate, ~ m^/0 2 ~ H, from which we deduce that 
Rd<t> = {m]p 15 <Po^ 5 Mp 15 /m^R^. At this value of R, we see that inflaton decays dominate the 
energy density only if O < m ^ 12 'Mp 3 /m^ 4 ~ 10~ 2 Mp. For larger values of the initial vev 
along the flat direction, reheating and thermalization is entirely determined by the dynamics 
of the flat direction [$]. In this case, we expect a baryon asymmetry of order 0(1) [TJ [2] 
unless it is diluted by other processes such as moduli decay or the vev is regulated by non- 
renormalizable or GUT scale operators [Hj. 

For relatively large </>o, the oscillations of the flat direction come to dominate the en- 
ergy density. In this case the decay of is determined by = (m ! j ) /(j) 2 ) m 2 b )(R/R t p) 3 = 



(m 1 p(j) /Mp)(R i p/R) 3 ^ 2 = H. The scale factor at decay is R^ = (<pQm^ 3 /m 1 ^' 9 Mp)R, 



Because this decay occurs relatively late, the reheat temperature is suppressed, Tp ~ 



The above argument for the decay of ignores the possibility that the flat directions decay 
on a much quicker timescale in a nonperturbative fashion (preheating). In the following, we 
show that these effects can be indeed dominant, and lead to the quick depletion of the flat 
directions within their first (~ 5 — 10) rotations. If so, the Universe remains dominated 
by inflaton oscillations until inflaton decay. Reheating in this case is produced by inflaton 
decays, and the baryon asymmetry is expected to be large unless O is regulated by non- 
renormalizable operators or dilution due to moduli decay occurs. 

The plan of the paper is as follows. In Section [2j we discuss an argument, based on a toy 
model of preheating, which has been used to claim that preheating along flat directions is 
negligible. In the remainder of that Section, and in Section [31 we explain why this conclusion 

3 As we show below, nonperturbative effects associated with mixing dramatically alter this conclusion. 




(3) 





10 6 GeV. 
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might not remain valid when realistic models are considered. In Section HI we present an 
explicit computation of the actual nonperturbative decay. Section [5] contains our conclusions. 

2 Nonperturbative decay of the flat direction 

In this section, we will first consider the case of two complex fields with a global U(l) 
symmetry. We will see later that the time evolution of the mass matrix along the flat 
direction leads to preheating. However, when the U(l) is gauged, the fields involved in the 
preheating are simply the massive Higgs and its related Goldstone boson eaten by the gauge 
field. We next consider a case with two flat directions involving four complex fields. In 
this case, we find that in addition to the four real degrees of freedom associated with the 
decoupled flat directions, and the massive Higgs and Goldstones fields, we are left with two 
real degrees of freedom for which preheating is found to occur. 

We saw in the previous section that when preheating is neglected, the flat direction comes 
to dominate over the inflaton decay products, if its initial amplitude is sufficiently large, 
4>o > 10~ 2 M P . In this case, the thermal bath which is relevant for the following history of 
the Universe is the one formed at the decay of the flat direction. Most existing studies focus 
only on the perturbative (single quanta) decay of . The possibility of a nonperturbative 
decay along flat directions was considered in [T6l [T7] where an important difference between 
this case and nonperturbative inflaton decay at preheating was pointed out. In the latter 
case, the inflaton is typically assumed to be a scalar field oscillating about the minimum 
of its potential, where its amplitude vanishes. Flat directions are instead complex fields. 
Nonrenormalizable operators (as well as susy breaking masses), give rise to terms in the 
potential which also depend on the the phase a of the flat direction ((f) = \(f)\ exp (ia) ). This 
in general leads to a time evolution of o , which can be seen as out of phase oscillations 
of the real an imaginary parts of the flat direction; hence, <fi never vanishes during these 
oscillations. 

The general picture was originally discussed in [TJ, where the "angular" motion of <p 
is associated to a baryon number charge. A nonrenormalizable operator is responsible for 
an initial a ^ . Subsequently, the field evolves to a lower amplitude, where the dominant 
potential term m^\(f)\ 2 is a independent. This results in an approximately elliptical motion of 
the flat direction in its complex plane (the amplitude slowly decreases due to the expansion 
of the Universe, so that <p is actually spiraling down towards the origin). For the following 
estimates, we can assume that the amplitude of <fi oscillates between the maximum O and 
the minimal value e0o over a timescale mT 1 . The value e = corresponds to straight line 
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radial motion of <fi in the complex plane (a = 0); in this case the phase a can be rotated 
away, and 4> oscillates as a real field. In the supersymmetric case, the value of e depends 
on the actual potential, and on the initial conditions for <fi. Values 10~ 3 < e < 10 _1 are 
typically found [17] . 

In analogy to what has been computed for inflationary preheating, one can consider [161 
[T7] the excitation of a single (complex) field Xi due to its coupling with the background flat 
direction, 0, 

= g 2 \ct>\ 2 \ X \ 2 (4) 

This interaction leads to a time dependent effective mass for the real and imaginary compo- 
nents of x, 

m 2 eS , x = m 2 x + g 2 \^(t)\ 2 (5) 

(the "bare" mass m x is expected to be of order the electroweak scale, and can be neglected for 
these considerations). The nonperturbative decay of the flat direction takes place whenever 
the frequency of the quanta of \ varies nonadiabatically, Co > uj 2 . For relatively small 
momenta p, one finds to 2 = p 2 + m 2 s ~ g 2 \<fi (t) | 2 . The ratio Co /uj 2 is maximized when the 
amplitude of <fi is minimal. When this occurs, one finds Co/io 2 ~ ge<f)om r f ) /g 2 e 2 <j)Q . Therefore, 
4> decays nonperturbatively only if 

^ (6) 

Since, is of order the electroweak scale, while 0o is close to the Planck scale, even a 
very small deviation from radial motion invalidates this inequality. Based on this considera- 
tion [16j[T7], it was concluded [9] that the flat direction decays only perturbatively, leading 
to a suppressed reheating temperature and a solution to the gravitino problem. 

However, this conclusion strongly depends on the coupling (j4j) assumed in these analy- 
ses, and in the resulting mass term ([5]). In concrete cases, the coupling is typically more 
complicated. For instance, the interaction may couple the real and imaginary components 
of x j or could couple \ to some other field. Furthermore, the coefficients of these couplings 
may depend also on the phase a of the flat direction, in addition to its amplitude. As we 
argue in the next Section, this is actually what happens for realistic MSSM flat directions. 
For definiteness, consider a flat direction involving two complex scalar fields with a global 
U(l) symmetry and a potential 

V = 9 j (\$i\ 2 - \®2\ 2 ) 2 (7) 
which has the same structure as the MSSM potential from the D— terms which we will 
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discuss below. If we expand the fields about their background value, 

$1 = + (f + X ) 

and 



we can write 



which when inserted into (171) becomes 



V = AV 



fr 



+ (Z-x) 



(4Re(0 X *)+4Re(rx)) 2 



(9) 



(10) 



This interaction couples the real and the imaginary parts of x (which we denote by x = 
(%i + 2X2) /v2)) and, consequently, the effective mass (jSJ) is now replaced by a 2 x 2 nondi- 
agonal mass term 

• 2 r ~ " ~ f 2 cos a sin a \ 



= (xi , X2) t Mfi (Xi , X 2 ) 



A^ 2 = 



m; + f cos 2 a 
/ 2 cos a sin a 



where / = y/2g\(f)\. Neglecting the small electroweak masses m Xl , m X2 . 



one linear combina- 
tion of xii X2 is massless, while the other one has mass f 2 . Since this eigenmass only depends 
on the amplitude of <fr > one may be tempted to conclude that the above considerations still 
apply, and that nonperturbative effects are negligible in this well. 

However, the nature of the massive and massless combinations is time dependent; when 
a = the massive combination coincides with xi > while Xi is massless. The real and 
imaginary components interchange their role every quarter of rotation. We will see that 
this time variation is also a source of nonperturbative particle production, which, being 
proportional to &, is not suppressed as e increases. In particular, it takes place even if 
the eigenmasses are constant or vary only adiabatically. As we show in Section 14.21 strong 
nonperturbative production takes place even if the motion of <fi is perfectly circular, and \(f)\ 
is constant. 

Now if we gauge the U(l), as the forthcoming analogies with the MSSM will require, we 
can show that the massive combination of xii X2 is in fact a Higgs boson while the other 
combination is the Goldstone boson. To see this, consider the coupling between the gauge 
boson and the scalar field excitations emerging from the covariant derivatives of $j. This 
coupling could be removed by an infinitesimal U (1) gauge transformation, with the gauge 
parameter A satisfying 

r sin<7 xi — co s cr X2 



d. 



do 



x/2|0| 
sin a Xi — cos a X2 



y/2, 



a 



(cosctxi + sina^) 



(12) 
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In the standard higgs mechanism, a = 0, and one can completely decouple the gauge boson 
from the neutral excitations by choosing 

_ sin axi- cos a X2 

The unitary gauge eliminates the massless Goldstone boson and the remaining massive higgs 
excitation is decoupled. In the present case, f|T3|) is clearly not a solution of f|T2|) ; however, it is 
still convenient to make this gauge choice, since it eliminates the massless neutral excitation, 
leading to the following quadractic action for the fluctuations 

\D,H 2 ~ V ^\ ( 9 »Xm) 2 + \ (<V) 2 x 2 m + 2g 2 \<P\ 2 Al - Ag\ct>\aA oXm - 2g 2 \<\>\\ 2 m (14) 

where Xm = (cosaxi + sin 17x2) is the massive higgs, while the orthogonal massless scalar 
has been absorbed in the longitudinal part of A^. The mass matrix for the {<f) m , A Q } system 
is still nondiagonal, due to the time evolution of o . However, the nature of the corresponding 
eigenstates does not rotate in field space, and preheating effects can be neglected as long as 
a< |0| . 

Even in the gauged case, significant preheating from the rotation of the mass matrix will 
take place in presence of additional coupled degrees of freedom beyond the Higgs, Goldstone 
pair. This will happen in presence of 2 or more flat directions. To see this, consider a 
second toy model model with 2 flat directions characterized by four complex fields, $j with 
a potential 

V = 9 - (q |$i| 2 - q |$ 2 | 2 + q' |$ 3 | 2 - q |$ 4 | 2 ) 2 (15) 

where we have assigned the fields $1^ and $3,4 charges ±q and ±q', respectively. We expand 
$1,2 as above, and 

$ 3 = 4>' + (£' + x') , $4 = 0'+(£ / -x / ) (16) 

and, for simplicity, we take q'\(f>'\ = q\4>\ ■ In this case, the part of the potential quadratic in 
the fluctuations reads 

P 2 
/W = — (cos (7Xi + cos °' X2 + sin o- Xi + sin a' x' 2 ) (17) 

(where a' is the phase of </>'). 

The corresponding mass matrix has one massive and three massless eigenstates, whose 
nature changes with time in the {xi, X2, Xii x'2) basis, as the two flat directions rotate. We 
have verified explicitly that all four states are produced at preheating (the computation is 
analogous to the one for the mass matrix (jlip . which is described in Section 4), as long as 
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a ^ ±cr' . When the U (1) symmetry is gauged, one massless combination (the Goldstone 
boson) is eliminated; however, the three remaining eigenstates still rotate in field space, and 
the corresponding eigenvalues are produced during preheating. In the special case when 
a' = a (a 1 = —&), we are effectively back to a single flat direction, and indeed we note that 
only the two combinations \i + Xi an d X2 + X2 ( X2 — x'2) are entering in AV; in this case, 
the result (or lack thereof) for preheating is identical to the one for the toy model (J7J). 

We conclude this Section with a note on the nonlinear interactions of the scalar modes 
and on their possible effects on preheating. The excitation of fields produced at preheating is 
typically encoded in a large variance for these fields. However, if for example a self-interaction 
A% 4 is present in the potential, then energy considerations preclude the formation of a large 
variance. Therefore, we would conclude that such a field % cannot be significantly excited 
during preheating. This is particularly true if A is as large as typical gauge couplings [19|. 
These considerations, however, do not affect the toy models just discussed. From eq. (fTOl) 
we see that the nonlinear interactions of x involve the field £ , which is not produced at 
preheating. Energy considerations then only limit the growth of the variance of £ (which 
could be produced when rescattering effects of x are taken into account), but (x 2 ) itself can 
have significant growth, as long as £ remains sufficiently small. An analogous situation takes 
place in the second toy model, where the nonlinear interactions involve the fields £ and £' 
which are not produced at preheating. 

3 The mass matrix for fluctuations about a flat direc- 
tion 

The mass matrix M 2 introduced in the previous Section (namely the one in eq. (1111) without 
the small susy breaking masses m X i) is precisely the one found for fluctuations around a 
MSSM flat direction, when one considers the part of the potential coming from D— terms: 

a 2 

V D = —D a D a , D a = X*T a Xi (18) 

where the index a runs over the generators of the gauge group, and Xi are the MSSM fields 
charged under that group element. To find the mass term, consider the MSSM excitations 
around a generic flat direction with vev . In this way, each MSSM field Xi is either of 
the form <p + 5Xi (if the field takes part in the flat direction) or (DQ. By definition, each 
D a vanishes when evaluated on a flat direction, that is when the fluctuations 5Xi are set to 
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zero. For this reason, D a is at least linear in the excitations, 
D a = <f>* J2 CiSXi + h.c. + O (5X 2 ) 



= 2\4>\ 



coscr 



Re (V^ Ci 5xA + sin a Im ^ aSxA + O (dX 2 ) 



(19) 



where c, are numerical coefficients dependent on the flat direction and on the generator 
considered^. It is clear that squaring (fT9l results in general in a mass matrix which depends 
both on the amplitude and the phase of 0, and which is non-diagonal in field space. Moreover, 
if we go to a basis constructed from the linear combinations entering in the linear term, the 
resulting mass matrix will be block diagonal, with all nonvanishing blocks proportional to 
M 2 . 

3.1 A single flat direction 

It is instructive to see this explicitly for a specific example. We consider the LLE C flat 
direction, and for definitness, we give a vev to 



Id 



= Iji) = (t c ) = <f (20) 

In addition, we define 5X = (5Xr + i5Xj) / \/2 to be the fluctuation of each MSSM fie 
about this flat direction. The mass matrix from the D— terms involves 10 coupled real field 
They are the real and imaginary parts of the fluctuations of the three fields involved in the 
flat direction (system /), plus the real and imaginary parts of 5e and (system II). The 
two systems are actually decoupled at the quadratic level, since the mass matrix for the fields 
in the first system comes from the diagonal D 3 (isospin) and D' (hypercharge) D— terms, 
while the one for the second system arises from the SU(2) off diagonal generators. 
The off diagonal structure for the first system arises in the 

f fa ~ V)r (^e ~ A*) I (^ e + ^-2T% (v e + fj, - 2t c ) j (v e + fj, + T c ) R (v e + fj, + T°) f \ 

\ V2 ' V2 ' V2 V2 V2 ' y/2 J 

(21) 

basis, where the mass matrix reads 

g 2 \<j)\ 2 M 2 
M 2 = | 3g' 2 \<j)\ 2 M 2 j (22) 




The linear term (|19p arises from mixed terms in (| 18[) . in which we take the background <j) from one of 
the two Xi and the excitation from the other. 

5 We are focusing here only on the part of the potential coming from D— terms. The part of the potential 
from the F— terms gives a mass matrix for a separate set of MSSM excitations, and it is diagonal in field 
space and cr— independent. 
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We note that each entry in this mass matrix is actually a 2 x 2 matrix. The first nonvanishing 
block is the mass matrix for the coupled system composed by a heavy higgs and the Gold- 
stone boson which provides the longitudinal component of W3. The second block is instead 
the mass matrix mixing another heavy higgs with the goldstone boson which provides the 
longitudinal component to B. The remaining two massless eigenstates are the actual real 
and imaginary components of the flat directions, and are decoupled at the quadratic level. 
For the second set of fields, in the basis, 



where each block is the mass matrix between a heavy charged field and the Goldstone boson 
providing the longitudinal component to W± and W2 , respectively. 

In summary, the LLE C flat direction breaks SU(2)x U(l) completely. For each broken 
symmetry, there is a massive (higgs) and a massless (Goldstone) excitation; in addition, the 
are two massless excitations (whose nature does not change with time) which represents the 
actual fields associated with the flat direction. This pattern (1 massive and 1 massive field 
for each broken symmetry + 2 additional massless fields) applies for each flat direction. 

This shows that preheating from mixing does not occur if we have a single flat direction. 
The field redefinitions (12T1) and (1231) are constant in time, and do not lead to particle produc- 
tion. After these redefinitions, we are left with a sum of systems analogous to the first toy 
model discussed in the previous Section, where the quick rotation in field space is eliminated 
once each goldstone is absorbed in the longitudinal component of the corresponding gauge 
boson (cf. eq ( TH1) .) However, the quick rotation of the mass matrix in field space, and the 
corresponding preheating effect, takes place if two or more flat directions are excited, as we 
have discussed at the end of the previous Section (cf. the discussion after eq. ffTTl) ). 

3.2 Two or more flat directions. 

The flat direction considered above allows several other independent flat directions to be 
excited simultaneously. Note that this is not a general feature of all flat directions. For 
example, the H u Hd flat direction makes it very difficult to excite other directions which 
remain F-flat. However, there are many mutually non-exclusive directions, including the one 
discussed in the previous subsection. 




(23) 



we find instead 




(24) 
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For example, let us consider the simultaneous presence of the two flat directions LLE C 
and QLD C . For definitness, we give vevs to 

Ou) = |0l (r c ) = \<P\e td (ik) = W\ (sl) = W\e ia {ve) = JW + W\ 2 (25) 

One can check explicitely that the scalar potential is both F- and D- flat along this direction. 
Of the 5 potential phases, only the combinations (arg^ e + arg/x + argr c ) and (argz/ e + 
argfii + args^;) are gauge invariant. Thus we are free to make the above assignments which 
will prove useful below. This combination of vevs will break the standard model gauge 
group SU(3) C x SU{2) L x Z7(l)y — > SU(2) C , ie. leaving a color SU(2) subgroup unbroken. 
As a result, we expect 9 massive states corresponding to Higgses, and 9 Goldstone modes. 
Our system is clearly more complicated now. Nevertheless, in analogy with our previous 
discussion, we can separate the excitations which participate in the mass matrix as follows. 
The real and imaginary parts of the pairs (o?2, Sj) and (o?2, S2) will split into two 4x4 blocks 
as in eqs. ( 1231 and fl24l) . each containing two massive Higgses and two Goldstones responsible 
for the longitudinal components of the four off-diagonal color generators (Ti j2 and T 4j5 ). 

Next, we can further isolate a 6 x 6 matrix for the fields, e#, ej, Vn R , Vm^i, ui R , Mi/.This 
matrix contains 2 Higgses which break the off diagonal generators of SU{2)i J and their 
corresponding Goldstone modes. Due to our choice of phase assignments, the entire matrix 
is phase independent and the two remaining massless fields do not contribute to preheating. 

The remaining 10 fields, are the fluctuations of the fields with nonvanishing vevs, namely 
v eRi Vei-, Hr, Hi, t c r, t c /, di R , d u , sl R , s\j. Their mass matrix arises from the D— terms cor- 
responding to the broken gauge symmetries with diagonal generators (namely, hypercharge, 
isospin, and the sum of T3 and T§ color generators). For this reason, they form a system 
which is analogous to system / described in the previous Subsection. Specifically, the three 
states v ei , fii, and d u do not enter in the mass matrix, due to the above phase assignment for 
the flat directions. The mass matrix for the remaining 7 fields depends on the phases a and 
9. The overall system of 10 states is characterized by 3 massive eigenstates, 3 Goldstones 
(which provide the longitudinal components to the gauge fields of the broken diagonal sym- 
metries), and 4 light physical states. Clearly, all the eigenstates are linear combinations of 
the 10 states of the system. What is relevant for the present discussion, is that each of the 4 
light physical states has some coefficients in this linear combination which are proportional 
to the cosine and sine of the two phases a and 9 . Therefore, these states are produced by 
preheating from mixing, in an analogous way to what happens in the two flat direction toy 
model of the previous Section. 

It is interesting to note that the number of fields excited may be larger still. For example, 
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a direction involving U C D C D C with (c^) = (dg) = (b\) = <p" is compatible with the previous 
two flat directions discussed. While the remaining SU(2) C subgroup is now broken, The 
new direction involves nine new complex fields. Of these 18 new degrees of freedom, three 
are Higgses, three are Goldstones, two are decoupled, but 10 are light and are involved in 
mixing. 

We note that mutually non-exclusive flat directions are not restricted to simple monomials 
involving three chiral superfields. The direction LLD C D C D C is compatible with QQQL. For 
example, giving vevs to u e , d\, s%, 63, u 2 , d 3 , Ci, and r would involve 54 fields in the mass 
matrix. With 12 Higgs and Goldstone modes, we are left with 30 fields which may participate 
in preheating! 

We also wish to stress that these flat directions are expected to be excited by scalar 
field fluctuations during inflation. In that case, all mutually non-exclusive directions will be 
excited. In other words, if the LLE C direction we discussed in the previous subsection is 
excited, there is no way to prevent the other compatible flat directions from acquiring large 
vevs as well. 

As we did at the end of Section [2j we conclude also this Section with a note on nonlinear 
interactions. Contrary to what happens in the toy model with two flat directions, we expect 
that, for a generic flat direction, the fields produced at preheating have large (namely, of 
gauge strength) self-interactions, as well as large interactions with some of the other MSSM 
fields. It is well known from models of inflaton preheating that large self-interactions can 
strongly limit the excitation of a field at preheating. For instance, if we couple a field \ to 
the inflaton ip through a quartic g 2 ijj 2 x 2 interaction, strong preheating is expected through 
parametric resonance. However, the production of x is strongly suppressed in the presence 
of a quartic interaction Ax 4 , with large coupling A. Unfortunately, the self interactions 
around realistic flat directions are much more complicated. Each D— term contributes a 
term of the form D a D a ~ (vev x linear + quadratic) 2 to the potential, where linear and 
quadratic refer to a sum of terms of that degree in the fluctuations (cf. eq, ( 1T91) for a more 
detailed expression). It is plausible that variances of different fields can grow significantly 
at preheating, but that the various terms entering in this expression, once summed, will not 
contribute significantly to the higher terms (more than quadratic) in the potential. Moreover, 
the presence of cubic terms in the potential can lead to two possible effects. Firstly, some 
MSSM fields may also develop a large vev, as the variances of some other fields are growing 
at preheating. Secondly, these terms lead to interactions which correspond to the scattering 
of a quantum against the zero mode and could enhance the decay of the flat direction. These 
considerations can strongly affect the nonperturbative decay of the flat directions, and they 
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can presumably be studied with the aid of numerical lattice simulations [18J. 

4 Computation of the decay 

In this Section we return to our original toy model and we show, by explicit computation, 
that complex flat directions can decay very quickly during their first few rotations. We 
have shown that in more realistic cases, the mass matrices take on a very similar form. We 
divide the presentation in three subsections. In the first, we present the general formalism 
needed for computing the decay. We then perform a numerical computation for a particularly 
simple background evolution, chosen to single out the specific effect we are considering. In 
the third subsection, we discuss the computation in a more realistic context. There we will 
also comment on the case toy model with two flat directions as we have seen this to be 
required by simple counting of degrees of freedom. 

4.1 General formalism 

To compute the nonperturbative decay of the flat direction 0, one first quantizes the fields 
Xi coupled to it. The quantization occurs in a time- dependent background, due both to the 
expansion of the Universe - encoded in the scale factor R - and to the evolution of <fi (treated 
as a homogeneous background field). The time evolution generates quanta of the fields Xi > 
at the expense of the energy density stored in the coherent motion of <ft ■ The standard 
computation is typically performed with only one field, x- When more fields are present, 
there are additional terms, which, as we see below, can be the source of very quick particle 
production. The general case has been worked out in the second ref. of [20]; here we just 
present the final results, referring to Section 2 of that paper for a detailed and self-contained 
derivation. 

The system of rescaled fields Rxi has the (squared) frequency 

&ii (26) 

where k is the comoving momentum, prime denotes a derivative with respect to conformal 
time 1] (drj = dt/R) and £ controls the coupling of <p to the curvature (the value £ = 
corresponds to minimal coupling, while £ = 1/6 to a conformal coupling). In general, the 
frequency is not diagonal in field space, due to the nondiagonal mass term. However, we can 
rotate it through 

C T (77) n 2 (77) C (77) = u 2 (77) diagonal (27) 



^ = R 2 M 2 J + 



Till 
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We denote by \ = C T x the fields in the basis in which the frequency matrix is diagonal. We 
also denote by uf the i— th entry of the diagonal matrix uj 2 . The set of uji represents the 
frequencies of the (time dependent) physical eigenstates of the system Xi ■ 

Initially, Q 2 is constant or evolving adiabatically; the fields Xi are Fourier decomposed 
in terms of only positive energy terms, corresponding to a vanishing occupation number. 
The time evolution of the frequency also generates negative energy terms; the presence of 
both positive and negative energy terms signals the production of quanta Xi ■ 111 the single 
field case, the coefficients in front of the positive and negative energy terms, known as 
Bogolyubov coefficients, are denoted by a and (3 , respectively. For scalar fields, they satisfy 
\a\ 2 — \f3\ 2 = 1 . The occupation number of the the produced field is given by \j3\ 2 (which is 
initially vanishing). 

In the computation of [20], the Bogolyubov coefficients are promoted to matrices in field 
space. They are determined by the initial conditions a = 1 , (3 = 0, and by the evolution 
equations 

, u' 

a = —iua H (3 — la — J {3 

2uj 

(3' = —a + tuj(3-Ja-I(3 (28) 
2cl> 



where we have defined 



J, J= I l^T^=±^=T^A , T = C T C (29) 



2 \ 

(the upper and lower signs refer to I and J, respectively). Finally, the occupation number 
of the i— th eigenstate x% can be shown to be 

n l (r ] ) = ((3*(3 T ).. (30) 

It is instructive to comment on the evolution equations (|2"B1 . We observe that (3 remains 
zero as long as Q is constant. The terms proportional to u'/uj are instead the source of 
particle production in the single field case, where I = J = 0. They account for the well 
known gravitational particle production (due to the change of the eigenfrequencies u){ with 
the scale factor R), and for preheating (due to the change of the mass matrix M 2 (0) ). These 
terms can be neglected as long as w' < w 2 , which was the criterion used in the literature to 
claim that the nonperturbative decay of flat directions is negligible. However, in the multi- 
field case, there is an additional source of production, due to the variation of the original 
frequency matrix Q ; this leads to the terms proportional to /, J in eqs. (I28p . We show below 
that these terms can lead to the fast nonperturbative decay even if the eigenfrequencies u>i 
are constant or slowly varying. 



14 



4.2 Production in a simple background 

We now apply the formalism given in the previous subsection to a very simple background. 
We consider the interaction ( TTDT) between the flat direction and a (complex) field x-, which 
leads to the mass matrix (JTTJ). In the next subsection, we discuss the self-consistent and 
realistic evolution for in a FRW spacetime geometry. Here, instead, we treat them as 
an external background, chosen so as to single out the production of quanta of x from the 
rotation. Specifically, we work in a Minkowski spacetime, and we take a constant amplitude 
|0| ; in this way, the {xi, X2} system (real and imaginary components of x , respectively) has 
constant eigenfrequencies U 



ui = y/k* + f* > ^2 = k (31) 

and particle production is solely due to the varying phase of <fi (encoded in the matrices I, J 
in eq. ([28])). 

The mass matrix Ai 2 of the system is given in eq. (ITTj) . where a denotes the phase of . 
We denote by n\ and n-i the occupation numbers of the massive and massless eigenstates, 
respectively. According to the above discussion, these number densities have an unambiguous 
meaning only as long as Ai 2 is constant (or, at least, adiabatically evolving). For this 
reason, we assume that a varies significantly only for a finite interval of time, and that it 
slowly approaches two constant values at early and late times. A simple function with this 
property is 



o(t) = ixN 



1 + tanh 



(32) 



.irN. 

With this choice, the flat direction performs N complete rotations between t = — 00 and 
+00 . Moreover, a = at its maximum, which is precisely the value that it would assume 
if it was rotating with constant amplitude in the potential m^|0| 2 . We interpret rii (t — > 00) 
as the occupation number of the i — th eigenstate generated after N rotations of the flat 
direction. 

To find the occupation numbers, we performed a numerical evaluation of eqs. (|28l) . There 
are two mass scales in the problem given by io\ ~ / and by o' ~ . It is convenient to 
define 

/i^ = -^- (33) 

In realistic cases, /x can be as small as ~ 10~ 14 . The timescale for particle production is 
set by the evolution of <p , that is mT 1 . However, the intermediate matrices a, (3 evolve 



6 We assume here that the mass for the field x arises only from the interaction with the flat direction, that 
is rrixi = in eq. (JTTJ) . We discuss the effect of nonvanishing "bare" masses at the end of this subsection. 
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on the much quicker timescale f" 1 , dictated by the largest eigenfrequency. Therefore, the 
computation becomes progressively more time consuming as \i decreases. We performed 
evaluations for different values of /i , with values of /i as low as 10~ 5 . The results obtained 
exhibit a clear scaling with /i , so that the amount of produced quanta can be easily inferred 
beyond the range we have directly probed. 



1e+10 



100000 



1e-05 
1e-10 



-4 -2 2 4 

H^t/27i 

Figure 1: Growth of the occupation numbers for 2 rotations, for the specific choice [i = 10~ 3 
and p = yu/10. The occupation numbers are dimensionless since they represent number 
densities with respect to both coordinate and momentum space. 

We start the discussion of the numerical results by comparing the growth of the two 
occupation numbers for a given choice of parameters. We fix /i = 10 -3 and we consider 
modes with momentum k = m^/lO (for later convenience, we define the rescaled momentum 
p = k/ f). The evolution of the two occupation numbers n\z is shown in fig. HJ 

In the numerical computations we performed, the argument of the tanh in (|32|) ranges 
from —xq to +xq . The larger xq is, the more adiabatic the initial and final states are. The 
results shown in fig. [1] have been obtained for xq = 5 . The occupation number for the 
massless eigenstate grows during the entire evolution, and saturates at a constant value at 
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late times as seen in the figure. In contrast, the occupation number of the massive eigenstate 
is initially growing, but it eventually decreases as the flat direction comes to a rest. We found 
that Hi continues to decrease for larger values of xq. Since the mass of the massive eigenstate 
is much greater than that of the flat direction, we are tempted to argue that the production 
observed in the figure is an artifact of not having a perfectly adiabatic initial and final 
state, so that n\ should be vanishing in the limit of xq — > oo . This issue, although very 
interesting per se, does not affect our estimate of the decay time of the flat direction. For 
any finite value of Xq , the numerical results give an upper bound on the amount of massive 
quanta produced. We verified that even for x = 5 , the energy density produced in the 
massless quanta is several orders of magnitude greater then one obtained for the massive 
state. Therefore, we simply concentrate on the massless quanta, for which the result is 
reliable (since n 2 has saturated to a constant value). 
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Figure 2: Spectra of massless quanta produced after N rotations of the flat direction. 



In fig. El we show the spectra of massless quanta produced after N rotations of the 
flat direction. The spectra are characterized by an exponentially growing resonance band 
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of momenta lower than the mass of the flat direction, k < . The spectra shown in the 
figure have been obtained for /i = 10 -3 . We performed an identical study for several values 
of p , from 10 -2 down to 10~ 5 . We found that the spectra in the resonant band are only a 
function of iV and of the ratio pj p . Therefore, the energy density in the massless quanta, 
Px2 = I cPk k n 2 , scales as p 4 , and it grows exponentially with N . The energy density in 
the flat direction is instead ~ m 2 |0| 2 ~ p 2 f 4 / (2g 2 ) . From our numerical results, we find 



This result is reliable only as long as r x <C 1 , so that backreaction effects can be neglected. 
If backreaction is negligible until the end of the process, setting r x = 1 in fl34|) provides a 
good estimate for the decay time of the flat direction. This occurs after 



rotations. So, even for realistically small values of p , the nonperturbative decay takes place 
within a few rotations, on a much shorter timescale than the perturbative one, which would 
typically last ~ 10 11 rotations in a Universe dominated by the radiation from inflaton decays. 

As in the case of inflationary preheating, the most important backreaction effect is prob- 
ably the scattering of the produced quanta against the flat direction. This effect, commonly 
known as rescattering, destroys the coherent motion of 0, and - consequently - it switches 
off the resonant production of the quanta of \ . So, backreaction decreases the particle pro- 
duction below the estimate ([Mil . However, it also contributes to the depletion of the zero 
mode of 0, so that effectively eq. ( 1351) may still be a good estimate for the decay time of the 
(coherent) flat direction. We comment more on this issue (and on the following evolution of 
the system) in the concluding Section. 

Finally, we note that a nonvanishing "bare" (that is, <fi— independent) mass for \ can be 
readily included in this computation (we expect it to be present as a soft susy breaking term). 
A nonvanishing m Xl = m X2 = m x in eq. (TTTj) amounts to replacing k 2 with k 2 + m 2 in the 
computation. So, the parameter p appearing in the above figures is related to the physical 
momentum k as k = J (p f) 2 — m 2 . As a consequence, the spectra shown in fig. [2] should 



start at p — m x / f rather than at p — 0. The resonance band exists for y (p 2 + m 2 .) < rn^, 

or p < -y/(?ri 2 — m 2 ) = p. The resonance band disappears if m x > . When m x < , the 
resonance band shrinks to < p < p (with p < m$). While it is smaller quantitatively, we 
expect a qualitatively similar result, unless the masses are nearly degenerate (m x less than 
but almost equal to m^). 




(34) 




(35) 
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4.3 Production in a cosmological self-consistent background 



We now compute the production of particles in a self-consistent evolving background, char- 
acterized by the potential 



V 



+ ^ (0 4 + h.c.) + g 2 + h.c.) 2 + ^m 2 x \ X 



1 



(36) 



The flat direction starts to move when the Hubble parameter H decreases below its 
mass. The evolution is not purely radial towards the origin, but the flat direction acquires 
some angular motion due to the a— dependent term, proportional to A . In the Affleck-Dine 
mechanism for baryogenesis [1] , the quartic coupling of interest is found from a 1-loop baryon 
number violating interaction with two supersymmetry breaking mass insertions. As a result, 
we expect that A ~ g^m^/^l [U [21] 

In the computation, we rescale <ft — <Po F exp (ia) / R , where 0o is the initial amplitude of 
the flat direction (we set R = F = 1 as initial conditions). We use conformal time, in inverse 
units of / = \/2g<pQ , and, analogously to \i and p defined above, we define /i x = m x j f . In 
this notation, the flat direction and the scale factor R evolve according to 



F" + 



E>2, ( 2 _/2 JX 

K fj, - a - — 



F + AF 3 cos (4a) 







F' 



XF 2 sin (4a) = 



IT 



p/2 ( 



+ 



R 2 Pip 

PM 2 



(37) 



where is defined in eq. (PQ) and dominates the Hubble expansion so long as O & Mp. In 
(I37D, we have defined f p = <f) /M p , A = X/2g 2 . 

We evolved the background equations numerically, starting with at rest due to Hubble 
friction. It is convenient to define a rescaled Hubble parameter as H = f h , from which we 
define the parameter, JF: 

hi = TyL (38) 

where hi is the scaled Hubble parameter after inflation. Using the mass scales described in 
section 1, ~ 10~ 7 Mp and ~ 10 _16 Mp, we expect T ~ 10 9 . The larger T is, the 
longer is the initial stage, in which <p is frozen, and the field \ is only excited due to the 
gravitational expansion. Since our focus is on the productions of x due to the motion of the 
flat direction, in the following we present results for T = 10 2 , so that <ft can be consistently 
started at rest, but the gravitational production of \ is limited. However, we also discuss 
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below how the final results changes when T is increased. For the other parameters, we chose 
initial conditions leading to sizable rotation. Namely, for any /i , we chose 



A = /i 2 , f p = 0.1 , a in = 0.2 , T = W 2 (39) 

We remind the reader that it is possible that non-renormalizable operators lead to values 
of / (and fp) which are smaller than that assumed for the present calculation. In geneneral, 
smaller vevs will lead to a more rapid dissolution of the condensate as the effective value of 
H will be larger. Furthermore, as discussed in the introduction, for vevs significantly below 
the Planck scale, the flat directions never dominate the Universe even in the absence of 
non-perturbative effects. 

We can obtain approximate analytical solutions for the background, using the fact that 
(i) the inflaton oscillations lead to a matter dominated evolution of the scale factor, (ii) that 
the amplitude of <fi , decreases as i?~ 3//2 during the rotations, and (iii) that the time derivative 
of the phase is approximately equal to . This leads to 

TT2/3 

R ~ (3irFN) 2/3 , o' ~ fi (3nTN) 2/3 , F ~ m (40) 

(37riV) ' 

where N is the number of rotations of the flat direction in its complex plane (we will often 
use N as our "time variable", since this leads to a more immediate interpretation of the 
results). These approximate relations are in good agreement with the numerical results. 

The mass matrix for the real and imaginary component of \ is again of the form (fTTj) . 
with m Xl = m X2 = \i x f . After diagonalization, one finds the two eigenfrequencies 



Y = ]Jp 2 +r 2 4 + ^^- 1 )+ f2 > j = yp 2 +^ + f^-i) (41) 

while the matrices /, J appearing in (1281) acquire the form 

7 _ a' ui + uj 2 ( -1 \ J _ a' ui - u 2 ( -1 \ . . 

/ ~ 2 x ajt^ v 1 / ' / ~ 2 v^~2 v - 1 ; 

In the case of minimal coupling (£ = 0), the R"/R term in uu 2 dominates initially in 
the range of momenta we are interested in, leading to lo\ < . This leads to gravitational 
production of the second eigenstate. Our focus here is not on the gravitational production, 
but rather on the creation of quanta of x from the motion of the flat direction, which is 
unrelated to £ ; moreover, the formalism that we are using assumes that all the eigenfre- 
quencies are real. For these reasons, we assume \ to be conformally coupled (£ = 1/6). 
This can be also considered a conservative assumption, since in this way we avoid part of 
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the gravitational creation (there is still gravitational production related to the mass term 
-R 2 /i^, in (HIT) ), so that a conformal coupling results in less particle production relative to the 
minimally coupled case. 

We start the analysis by considering a single mode with p = fi = 10~ 4 . For definiteness, 
we choose \i x = fi/2 . In fig. [31 we show the evolution of the adiabaticity parameters w'Juf . 
We note that both uj'Jcuf <C 1, corresponding to an adiabatic evolution of the eigenfrequen- 
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Figure 3: Evolution of the adiabaticity parameters uj'Juf for the two eigenstates Xi an d Xi ■ 

cies. For this reason, one may be tempted to conclude that the nonperturbative production 
of quanta of x is negligible. 

This conclusion would be correct in the single particle case, or if the mass matrix was 
diagonal, so that eqs. ff28l) would appear without the matrices I and J . To verify this, we 
incorrectly set / = J = in these equations, and we show in fig. H] the evolution of the 
occupation numbers n^j of the two eigenstates of the system. Indeed, we see from the figure 
that setting I — J — results in negligible production. In this case, the only source of 
production would be from the variation of the eigenfrequencies, which, as indicated by fig. [3l 
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Figure 4: The incorrect evolution of the occupation numbers xi an d X2 ■ The matrices I 
and J are set by hand to zero. Notice there is negligible production, in agreement with the 
adiabaticity shown in fig. [3l 

evolve adiabatically. However, when we include the matrices / and J in a manner consistent 
with eq. Q42D . we find significantly more particle production, as shown in fig. [51 

Having verified the strong effect associated with the rotation of the mass matrix, we can 
now turn to the computation of the total particle production. We perform the computation 
for different values of \i , ranging from fi = 10~ 4 to fi = 10 -7 , studying the production which 
takes place within the first N = 5 rotations of . The spectra of produced quanta show a 
clear scaling with respect to the parameters of the model: 



where the function u (normalized to 1 at small k ), is given in fig. [61 
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Figure 5: The correct evolution of the occupation numbers Xi an d Xi ■ Particle production 
stems from the matrices I and J (now included consistently). 



Among the different scalings summarized in eq. (|43p . the proportionality to /i is the 
most accurate. This leads us to conclude that this result can be safely extrapolated to 
lower (and more realistic) values of \x . With regard to the scaling with the number of 
rotations N, we observe the presence of a resonance band (characterized by exponentially 
growing occupation numbers) at relatively low momenta. The scaling (|43p is actually only 
accurate for the resonance band; however, this is good enough for the computation of the 
total production, since the quanta produced at higher momenta are completely negligible. 

In analogy with the Minkowski case studied in the previous Section, we expect the reso- 
nance band to be present for physical momenta smaller than the mass of (ft , or (in rescaled 
units), p/R < fi . Due to the expansion of the Universe, increasingly more comoving mo- 
menta enter the resonance band as the scale factor increases. However, the modes which are 
in the resonant band from the onset of the oscillations of are the ones which experience 
the most production. When the flat direction starts evolving, R ~ T 2 ^ (this immediately 
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Figure 6: Function u entering in eq. (|43|) . The two plots have been obtained by numerically 
evaluating the particle production for \i = 10 -6 , T = 100 , and N = 5 and by scaling the 
result as in eq. 



follows from the definition of J 7 ). For T = 100, which is the value at which the argument 
of u is normalized, this amounts to R ~ 20 . At the end of the first rotation, we find in- 
stead (numerically) R ~ 100 . Therefore, for JF = 100 , we expect a strong decrease of the 
occupation numbers in this range of p/fi , as fig. ([6]) clearly confirms. 

Finally, these considerations also explain the scaling of eqs. with respect to T , since 
T 2 ^ is essentially the normalization of the scale factor at the beginning of the oscillations 
of 4> i and the results should clearly depend on the physical momenta, rather than on the 
comoving ones. There is actually another smaller effect of T which we have not included in 
eqs. (H2D- Besides the overall "shift" of the spectra that we have just discussed, we also found 
slightly more production as T increases. This is due to gravitational production before <fi 
starts evolving (the greater T is, the greater is the value of the Hubble parameter at the start 
of the numerical simulation). The result given in (J32D nas been obtained for T = 100, but - as 
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we mentioned - T can be expected to be as large as 10 9 in realistic cases. Therefore, eqs. (j4~3l 
should be considered as a conservative estimate of the amount of particle production. We 
also performed numerical computations with T = 10 9 , finding that the total energy density 
p x increases by a factor of ~ 3 . This however does not change significantly our estimate for 
the moment at which the flat direction decays (see eqs. fj4~9l and fl50l) below). 

After this discussion, we can now compute the energy density of quanta of x produced, 
and its ratio with the energy density stored in the coherent oscillations of <fi . We estimate the 
decay time as the moment at which this ratio becomes one (with all the cautionary remarks 
already made at the end of the previous subsection). 

In the resonance band, and in the time interval we are interested in, the eigenfrequen- 
cies (HU) can be estimated as 
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This, together with eq. f|43|) for the occupation numbers, leads to 
d 3 k {ijji ri\ + lu 2 n 2 ) 
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the energy density of the 2nd state dominates in the range of iV we are interested in; the 
integral over momenta is numerically found to be J dpp 2 u ~ 0.03 /i 3 JF 2 . So, 
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Using the above estimates for F and R . 
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giving the ratio 
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Notice that T cancels in the ratio, since its role is simply related to the value of the scale 
factor at the beginning of the rotations. For g 2 ~ 0.4 , the ratio equates to 1 for 



From this we conclude that the flat direction decays long before the inflaton field. 

Finally, we return to the our enhanced toy model with two flat directions. Here, even 
in the gauged case, we expect mixing between four degrees of freedom (the other four are 
decoupled and represent the two flat directions). We recall however, that only production 
of the two light degrees of freedom which are orthogonal to the Goldstone mode can be 
responsible for preheating. Indeed, when we generalize the analysis detailed above, and 
when the two flat directions are characteriized by two independent backgrounds, = |0|e Jcr 
and 0' = |0'|e* 9 , we find that the occupation numbers, n%,n2, and are all produced as 
in Fig. [5j Furthermore, we have checked explicitely that when the two flat directions are 
degenerate, and = 0', n 3 = n 4 = and we return to the case where only the Higgs and 
Goldstone modes are involved. Thus we require two independent flat directions to be excited 
in order for preheating to occur. As we saw in section 3, this feature is rather common in 
the MSSM. 

5 Summary and Conclusions 

The scalar potential of the MSSM is characterized by several flat directions. This flatness 
may be lifted by supersymmetry breaking terms, nonrenormalizable operators, and/or super- 
gravity corrections. The latter are of particular relevance in the early Universe. For instance, 
if the MSSM fields have a minimal Kahler potential, then supergravity terms provide a mass 
(of order H) for the flat directions which prevents them from being excited during inflation 
[22]. However, if the Kahler potential exhibits a Heisenberg symmetry [13], such mass terms 
do not arise at the tree level and the flat directions can indeed develop a large vev during 
inflation [15]. No-scale supergravity [H] is one such example. The development of a large 
vev may have interesting consequences for the subsequent evolution of the Universe. 

Among the consequences of a large vev along a supersymmetric flat direction, is the 
possible generation of a baryon asymmetry [1] . Furthermore, the energy density stored along 
the flat direction may come to dominate over the inflaton decay products if the flat direction 
is long-lived. It is tempting to conclude that this is the case if one simply considers the 
perturbative (single quanta) decay of the flat direction. Indeed, the flat direction, due to its 
large vev, provides a large mass to the fields to which it is coupled, so that the perturbative 




(50) 
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decay proceeds by loop(s), with very heavy fields in the propagator. This typically results 
in a small decay rate. 

However, this consideration does not rule out the possibility that the flat directions 
decay nonperturbatively. It has been shown [23], that in many models of inflation, the 
inflaton decays nonperturbatively through a collective effect called preheating. In this paper, 
we have investigated the possibility of preheating during the evolution of fields associated 
with supersymmetric flat directions. In particular, it is important to determine whether 
preheating can lead to a significantly quicker decay than the perturbative estimates may 
suggest. Prior to our investigation, this question was addressed only in a very limited 
number of papers, which concluded that preheating is negligible in this case. Their claim, 
however, was based on toy models, where the flat direction <p is coupled to a single scalar 
X through a quartic |0| 2 |x| 2 interaction. The time varying amplitude of results in a time 
varying mass for the real and imaginary components of \ . However, in the specific case in 
which is a flat direction (rather than the inflaton) the variation of this mass is always too 
slow to lead to significant nonperturbative particle production. 

We have shown, however, that the couplings of the flat direction encoded in the MSSM 
potential are more complicated than the simple |0| 2 |x| 2 interaction. In particular, the part 
of the potential from the D— terms results in couplings of the form <p 2 XiXji leading a non- 
diagonal mass matrix for the fields x%- Even if the eigenmasses are slowly varying, the mass 
matrix M 2 itself varies with time. The dominant time dependence is due to the phase 
a of the flat direction, which is rotating quasi periodically (with a small decrease of the 
amplitude) around the minimum of its potential. The quasi periodic variation of M 2 leads 
to a very strong parametric resonance effect. We found that when at least two flat directions 
are excited, the flat directions typically loses most of their energy density within its first 
5 — 10 rotations due to these couplings. From this point on, backreaction effects of the Xi 
fields can be expected to play a very relevant role in the evolution of the system. Their 
complete investigation requires the use of numerical simulations on the lattice [18] . which 
are beyond the aim of the present workB However, due to the strong interactions between 
the MSSM quanta, these interactions are likely to lead to a fast depletion of the (remaining) 
vev of the flat direction. As such, it is very unlikely that the scalar field oscillations of 
the flat direction can ever come to dominate over the inflaton field and its decays product, 
particularly in the case in which the inflaton itself decays only gravitationally. 

Although the final stage of the decay remains to be explored, we have thus showed 

7 Numerical simulations are also required to determine if non-linear interactions affect the early stages of 
preheating. 
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that flat directions can potentially have a very fast decay channel which has been so far 
overlooked in the literature. There are several important consequences of this rapid decay. 
If flat directions never dominate the energy density of the Universe, entropy production 
and reheating are determined by the dynamics of the inflaton. The baryon asymmetry 
will be measured with respect to the entropy density produced by inflaton decays and will 
typically be large requiring either further dilution of the baryon asymmetry by later entropy 
production (e.g. due to moduli decays) or regulating the initial vev along the flat direction 
with some non-renormalizable operator [8j. Finally, in the absence of a large vev and hence 
large particle masses, inflaton decay and the subsequent thermalization also proceeds rather 
quickly [T2] , leading to a relatively large reheat temperature after inflation which could have 
dire consequences on the gravitino problem. 

Note added : After the present work was completed, a new paper which discusses the 
decay of flat directions [21] came to our attention. The first part of [21] is devoted to the case 
of a single flat direction. The analysis, and the conclusions, coincide with those already given 
in the present work. We comment on two of the remarks made in the second part (devoted 
to the case of two or more flat directions), although neither of them has been supported by 
a concrete computation. The first remark claims that flat directions may have hierarchical 
vevs (a small ratio r between the two vevs), so that the system may be effectively reduced to 
the single flat direction case. In general, the amount of produced particles will be suppressed 
by the ratio r. However, for any finite r, the number densities of the non-goldstone massless 
fields still grow exponentially with the number N of rotations of the flat directions. For 
the toy model with two flat directions, we have verified numerically that most of the energy 
density in the flat directions decays after N ~ 10 when r = (as opposed to iV ~ 5 

when r = 1). The second remark is that nonrenormalizable terms in the potential, which 
depend also on the phase of the flat directions, can in some cases result in chaotic motion for 
the phases, potentially suppressing the parametric amplification which takes place when the 
phases move with nearly constant velocity. First, nonrenormalizable interactions are very 
model dependent, and do not necessarily lead to chaotic motion for the phases. Second, 
even if they do, nonrenormalizable terms are relevant only for the first few oscillations of the 
flat directions. As soon as the flat directions evolve, their vev is decreased by the expansion 
of the universe, entering into a region where their potential can be safely approximated by 
the quadratic (phase-independent) part. When this happens, the phases will evolve with a 
nearly constant velocity, as the cases considered here. 
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